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A vacuum autoresonance accelerator scheme for electrons, which employs terahertz radiation 
and currently available magnetic fields, is suggested. Based on numerical simulations, parameter 
values, which could make the scheme experimentally feasible, are identified and discussed. 
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Lasers, pulsed; (140.7090) Ultrafast lasers; (350.4990) Particles. 


The aim of this Letter is to investigate conditions for an 
electron vacuum autoresonance accelerator scheme that 
would employ circularly polarized terahertz (THz) radia¬ 
tion (or T-rays) and currently available laboratory mag¬ 
netic fields, with the hope of stimulating future experi¬ 
ments. The system, subject of this Letter, is an electron 
(or electron bunch) injected in the common directions of 
radiation pulse propagation and an added uniform mag¬ 
netic field (see Fig. [Tj) . 

Autoresonance Laser Acceleration (ALA) of electrons 
has a relatively long history M- Theoretical vacuum 
ALA studies [3, 7] have shown that the electron would 
gain a tremendous amount of energy from the laser field 
when the initial injection energy, the laser frequency uj 
and the external magnetic field strength B s , all conspire 
to achieve resonance, or near-resonance. For axial injec¬ 
tion, and plane-wave fields, the resonance condition takes 
the form 


r _^ = eB ± /l+|o (1) 

UJ D TYIOJ V 1 - A, 

where, uj c = eB s /m is the cyclotron frequency of the elec¬ 
tron around the lines of the static magnetic field, and ujd 
is the Doppler-shifted frequency of the radiation field, as 
seen by the electron. Furthermore, m and e are the elec¬ 
tron’s mass and charge, respectively, and SI units have 
been used. Also, /3o is the speed with which the electron 
is injected, scaled by the speed of light c. Autoresonance 
essentially means the electron cyclotron frequency be¬ 
comes comparable to the Doppler-shifted frequency, seen 
by the electron, of the (circularly polarized) laser. When 
this happens, the velocity vector (3 of the electron and 
the electric field vector E of the laser field maintain the 
same angle during interaction, which leads, according to 
the equation 

ds 

- = - ec f3. E , (2) 

to synchronous energy gain by the particle from the ra¬ 
diation field. In Eq. m the electron’s relativistic energy 
is e = yrac 2 , where 7 = (1 — /3 2 ) -1 / 2 . On resonance, 
the vectors (3 and E gyrate, about the common direction 





FIG. 1. (Color online) Schematic of a possible ALA setup. A 
magnetic field B is used to bend the electron beam for collision 
with the THz pulse, which propagates along -j-z (propagation 
vector k). Front of the pulse catches up with the electron 
precisely at the origin of coordinates O. The magnetic field 
B s — B s k is responsible for the electron cyclotron motion. 


of the magnetic field and laser propagation, at the same 
frequency. According to Eq. the field-to-electron 

energy transfer rate is a maximum when (3 and E are 
antiparallel. 

To the best of our knowledge, a vacuum autores- 
onace laser accelerator has not been realized and none is 
currently available and in operation. Plane-wave-based 
studies have shown that, for the scheme to work, a strong 
magnetic field needs to be maintained over a long dis¬ 
tance, which makes such a device both expensive and 
prohibitively too long. Hope has been revived by recent 
investigations, in which petawatt optical-frequency laser 
pulses, modeled most realistically |8[ by Gaussian fields, 
have been employed. It has been shown that electrons 
can gain over 10 GeV of energy in a magnetic field of 
strength exceeding 50 T, and maintained constant over a 
distance of several meters Unfortunately, even these 
conditions are at present extremely difficult to realize. 

In this Letter, the laser fields will be replaced by THz 
fields, and the parameter space will be scanned in search 
of a parameter set which may be more likely to be exper¬ 
imentally realized in the near future. According to the 
condition m vacuum autoresonance may be achieved in a 
uniform magnetic field of a few tesla, provided the radia- 
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FIG. 2 . (Color online) Contour plot of the energy gained by 
a single electron, injected with 70 = 3 for interaction with a 
single circularly polarized Gaussian pulse, as functions of the 
magnetic held strength B s and the waist radius at focus wo- 
The pulse power is P = 100 TW, and its frequency is / = 4 
THz (A = 75 /im, period T 0 = 250 fs = FWHM). 

tion frequency is lowered by roughly two orders of magni¬ 
tude. Lowering the frequency by two orders of magnitude 
from the optical domain lands one in the THz region of 
the electromagnetic spectrum, roughly in the range 0.3 - 

10 THz. 

Low-power THz radiation has been the subject of in¬ 
tense investigation for many years now, and applications 
in such fields as imaging, wireless communication and re¬ 
mote sensing, have seen much progress recently [gHU. in 
pondering the idea of utilizing THz radiation in a vacuum 
ALA-like scheme, one must immediately come to grips 
with the need for high power. Currently considered the 
best source for high-power THz radiation, a free-electron 
laser (FEL) can gen erate radiation of only a few gigawatt 
(GW) power [10|,[l4|. Furthermore, promising candidates 
for generating more intense T-rays are the interactions of 
high intensity lasers with plasma or gas targets EM 2 . 
Thus, to experimentally realize a vacuum ALA-like THz 
scheme, the only remaining (and admittedly non-trivial) 
challenge would be to develop sources of THz radiation 
of power in the terawatt (TW) region and beyond, as will 
be demonstrated shortly. 

The single- and many-particle simulations in this Let¬ 
ter are all based on solving the relativistic energy- 
momentum transfer equations 

^ = — mP-E)-(E + c(3xB)}, (3) 

at 'jmc 

for the electron, subject to the adopted initial conditions 
of axial injection with scaled energy 70 = (1 — /3q) - 1 / 2 , 
together with the assumption that the front of the pulse 
catches up with the particle at t = 0 , precisely at the 
origin of coordinates. In Eq. ©, B is the sum of the 
radiation magnetic field and B s . Fields of the circularly 
polarized THz radiation are modeled by those of a short 
Gaussian pulse HMH. 



FIG. 3. (Color online) (a) Single electron energy gain (dashed 
blue) and energy gradient (solid red) as functions of its for¬ 
ward excursion distance. The parameters are those of Fig. [2j 
in addition to wo — 17A « 1.27 mm, and B s ~ 39.6 T. (b) 
Actual trajectory of the electron under the conditions of (a). 


In addition to the magnetic field strength B s and the 
THz frequency cj, a third parameter, namely /3q and, 
hence, the scaled injection energy of the electron 70 , plays 
a decisive role in determining the (exact, plane-wave- 
based) resonance condition (pQ). When a more realistic 
model is adopted, such as that of a Gaussian beam or 
pulse, detuning away from resonance results, and one is 
forced to search for near-resonance bv scanning the pa¬ 
rameter space of Fig. [2j for example jsj. 

For further discussion, we select from Fig. [2] a param¬ 
eter set which leads to near-autoresonance and, hence, 
high energy gain. With the gain defined by 

G = (7 - 7 o)«ic 2 , (4) 

Fig. [3pa) shows that an electron exits with kinetic energy 
Kexit ~ 396.2601 MeV from interaction with a single¬ 
cycle, 100 TW pulse, over a distance of less than 40 cm, 
and in the presence of a magnetic field of strength 39.6 
T. Note that most of the energy is gained over nearly the 
first 20 cm, and very little energy is gained beyond that 
point as the electron interacts with the weak tail of the 
pulse and is ultimately left behind. The energy gradient, 
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FIG. 4. (Color online) Exit kinetic energy distribution of a 
1000 ensemble of electrons. Interacting: Coulomb electron- 
electron interactions are turned on, and Non-interacting: 
Coulomb interactions are turned off. Electrons are assumed, 
initially, to be distributed randomly within a cylinder of ra¬ 
dius 0.232 /m i and height 4.642 gm, and centered at the ori¬ 
gin of coordinates. Electron kinetic energy is, initially, dis¬ 
tributed normally, with mean Ko = 1.022 MeV and spread 
AK 0 = 0.1%. 
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or the gain per unit forward excursion distance, 

dG _ f (3 - E\ 

dz~~ e \ f3 z )' 


(5) 


peaks over the first one centimeter and falls down to zero 
very quickly. The peak gradient is almost 8 GeV/m, 
or about 80 times the fundamental limit on the perfor¬ 
mance of conventional accelerators [20]. An average gra¬ 
dient , found simply by dividing the total gain by the to¬ 
tal forward excursion during interaction with the pulse, 
is about 1 GeV/m, in this example. 

Figure 03[b) shows the trajectory of the electron whose 
gain and gradient are shown in (a). The trajectory is the 
expected helix of increasing cross section Q. Note that 
the transverse dimensions of the helical trajectory are 
much smaller than the longitudinal electron excursion, 
which makes the trajectory essentially linear. 

The discussion so far has been limited to THz - ALA 
of a single electron. Next, we consider acceleration of 
an ensemble of 1000 electrons randomly distributed, ini¬ 
tially, within a cylinder centered at the origin of coor¬ 
dinates, and whose radius and height are 0.232 /im and 
4.642 fi m, respectively. Initial kinetic energy of the en¬ 
semble has a normal distribution of mean Ko = 1.022 
MeV and spread AKo = 0.1%. For the parameter set 
used in Fig. 03 our simulations yield the exit energy dis¬ 
tributions displayed in the histograms of Fig. 01 From 
the data, one gets K ex u = 396.2583 db 0.0028 MeV (or a 
spread of 0.0007 %), with the electron-electron Coulomb 
interactions turned off, and K ex u = 396.2562 zb 0.1682 
MeV (0.042 %), when those Coulomb interactions are 
properly taken into account. These many-particle results 
agree quite well with our single-particle calculations (Fig. 




FIG. 5. (Color online) Near-resonance magnetic field depen¬ 
dence (blue and circles) upon: (a) the THz pulse duration r 
in units of the radiation field cycle To, (b) the THz power, 
(c) the scaled injection energy 70, and (d) the THz frequency. 
For each data point in (a) - (d) the average energy gradient 
in GeV/m is given (red and crosses). Each data point is a 
result of calculations, along the lines of the work which led to 
Figs. [2] and 03 and employing their (other) parameters. For 
example, in (a) 70 = 3, the power is 100 TW, wo — 17 A, and 
f — 4 THz, and so on. 


03 K exit ~ 396.2601 MeV). Note also that the particle- 
particle interactions do result in a noticeable increase in 
the spread in exit energies. In absolute terms, however, 
the effect of incorporating the Coulomb interactions is 
small, due to the fact that the particle density of the 
initial ensemble is quite low ( r\j 1.273 x 10 21 m -3 ). 

Our focus in this Letter has been to demonstrate that 
THz - ALA of electrons may be experimentally feasible 
with present-day technology, as far as the needed mag¬ 
netic field strengths are concerned. However, the exam¬ 
ple presented above in some detail required the use of 
B s ~ 40 T. This kind of field strength is available only 
at large and expensive facilities [2l],[22j . For weaker fields 
to be utilized, other parameters have to be compromised. 
Values have to be picked for four more parameters which 
may ultimately lead to energy gains and energy gradients 
more modest than has been reported above. Figures [5][a) 
- (d) show how the values of B s , corresponding to near¬ 
resonance, vary with the pulse duration, the power of the 
THz radiation source, the scaled injection energy, and the 
THz frequency, respectively. Shown in red also are the 
corresponding average energy gradients, defined in each 
case as the total exit energy gain divided by the total 
forward excursion distance, during interaction with the 
pulse. The results displayed in Fig. [5] emerge from cal¬ 
culations following along the same lines leading to Figs. 
[2] and 03 for each set of parameters separately. In most 
cases considered, the average energy gradients are much 
higher than the natural limit on performance of the con¬ 
ventional accelerators, namely, about 100 MeV/m. Note 
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that in calculating these average gradients, the full excur¬ 
sion distance occurring during the full particle-field inter¬ 
action time, is used. Considering that most of the energy 
is gained during interaction with only a small fraction of 
the excursion, leads to the conclusion that the energy 
gradients are, effectively, much higher than is suggested 
by the averages reported in Fig. [5j 

Figure 0(a) shows clearly that the use of currently 
available laboratory magnetic fields in a THz - ALA 
setup is feasible, but the energy gradients that would 
be achieved fall to just a few times the performance limit 
of conventional accelerators Recall that the exam¬ 
ple considered above (Fig. [3]) made use of a one-cycle 
pulse. With several-cycle pulses, the resonance magnetic 
field strengths required fall already below 30 T, as shown 
in Fig. 0(a). The same conclusions may be drawn from 
Fig. 0(c). However, the opposite trend exhibited in Figs. 
[Sfb) and (d) confirms earlier conclusions. Increasing the 
source power and frequency calls for the need for stronger 
resonance magnetic fields, and leads to high gains and 


ultra-high energy gradients @,0. On-resonance depen¬ 
dence of B s upon the frequency / is linear in Fig. 0(d), 
essentially like the plane-wave case O. 

In conclusion, acceleration by cyclotron autoresonance 
of electrons in the simultaneous presence of THz fields 
and a uniform magnetic field of strength currently avail¬ 
able for laboratory experiments, is feasible, provided a 
THz radiation source of TW power is also available. This 
is the case when the electrons are injected along the com¬ 
mon directions of magnetic field and radiation propaga¬ 
tion. According to Fig. 0(c), even acceleration from rest 
is possible, but the magnetic field strength needed for 
such a scheme is about 300 T, in this case. 
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